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1. Introduction 

Effective mass theory had been used for years in several branches of modern physics like nuclear 
physics 01 or solid-state physics (see e.g. the pioneering works |2|-|4|). This theory is an useful 
tool for studying the motion of carriers in pure crystals and also for the virtual-crystal approximation 
to the treatment of homogeneous alloys (where the actual one-electron potential is approximated by a 
periodic potential), as well as of graded mixed semiconductors (where virtual-crystal potential is not 
periodic). The salient feature of this theory is that it approximates a complicated physical situation to 
the solution of a Schrodinger equation with position-dependent effective mass function, the so-called 
effective mass (EM) eigenvalue equation. The position-dependent EM is also used in the construction 
of pseudo-potentials, which have a significant computational advantage in quantum Monte Carlo 
method |5 1. Needless to emphasize that the growing interest in semiconductor physics parallel to the 
modern development in fabricating nanostructure technology creates a renewed attention to study the 
behavior of one-dimensional EM eigenvalue equation 

H EM {x)1>(x) = [T EM (x) + V(x)Uix) = Ei/>(x) (1.1) 

from a theoretical standpoint. In this equation Tem{%) is the EM kinetic energy operator and tp(x) is 
the EM wave function (also known as 'envelope wave function' in the literature |4|-|18|). It has to be 
mentioned that the justification of this EM approximation in the context of a realistic situation is not 
the purpose of the present study, rather we will concentrate on the analysis of the spectral properties of 
the EM eigenvalue equation The first step is certainly to choose a suitable form of the Hermitian 
kinetic energy operator arising from non-commutativity of momentum operator p = —ihd x and the 
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effective mass operator m(x). Different forms had been proposed in the literature, most of which 
may be written as a special class of the general two-parameter family proposed in Ref. 1 9 1 



with the constraint a + /? + 7 = —1 over the ordering parameters. 

Considerable efforts were made to remove the non-uniqueness of the kinetic energy operator 
II. 2\ or, in other words, to fix the values of the ordering parameters a, (3, 7. In |19| a step 
potential and a step mass were considered and it was shown that a = 7 is the only physical 
choice for an abrupt heterojunction. Later in an attempt to fix j3, two different conclusions were 
drawn, namely (3 = for a one-dimensional model |20| and (3 = — 1 for a three-dimensional 
model [21 1. On the other hand, in a series of works 1141 1221 l23l . the authors concluded that 
(1 = 7 = 0, (3 = —1 for an abrupt heterojunction. Among these works, Ref. |23| deserves 
to be mentioned separately because it presented the first example of a continuous function m{x) 
across the heterojunction. A new kind of kinetic energy operator was proposed 1151 for strained 
heterostructure, which is not included in (II .21 . in general, for position-dependent lattice constant. 
It should be mentioned that the choice a = 7 = 0, [3 = — 1 gives rise to the kinetic energy 
operator Tem = p(^/2m)p, which was first proposed in Ref. 0. Choosing the same operator 
some interesting pedagogical models were considered 11171 1241 to show the qualitative differences 
in quantum mecanical observables (e.g. reflection and transmission coefficient, band-structure, etc.) 
between EM and constant-mass case. Many other forms of kinetic energy operator had also been 
proposed, e.g., a = —1, (3 = 0, and 7 = which gives from (II. 2> Tem = (l/4m)p 2 + p 2 (l/4m) 
l25l : a = 7 = -l/2,/3 = which yields T EM = {l/2^m)p 2 {l/2^/m) |26|, etc. A different 
variation was derived, via path-integral formalism |27|, which comes from ll.2> for the following 
values: a = (— y/2 + i)/3\/2, (3 = —1/3, 7 = a*. It is therefore clear that no universal choice for 
the ordering parameters exists in the literature of EM theory. 

In recent times several authors [28-50] either started from a preferred ad-hoc choice for a, (3, 7 
or kept them arbitrary. In both cases attention had been paid to solvability for various smooth 
functional forms of V(x) and m(x) by employing the existing tools like supersymmetry [29, 
30, 34, 35, 40, 44-^8], Lie-algebraic approach 02] [33] |35l |39), shape-invariance |43] |46|, etc. 
The connection between solvability and the ordering parameters in equations J 1 . 1 1 and ll.2l > was 
discussed in Ref. 1311 . On the full line, smooth functions (in the sense that ml and m" are also 
continuous) were chosen for the first time in |28|, where the authors however concluded again that 
a = 7 = 0,/3 = — lby comparing their results with a limiting case where the potential and mass 
become abrupt. It may be mentioned that in several works [33-35] mass function was kept arbitrary 
and thus the solutions provided there were only formal. On the contrary, in most cases where smooth 
functional forms were chosen for m(x), we notice that m(x) — > as \x\ — ► 00. One possible way of 
eliminating this nonphysical situation is to consider the variation of the mass in a finite region (being 
constant outside), which is also natural on realistic ground. Of course, this correction will force us to 
obtain correct matching conditions for the model. 

Quite justifiably we will use the generic kinetic energy operator II. 21 without imposing any 
additional constraint over the ordering parameters. One of the purpose of our present work is to 
show that a sensible quantum situation could be modelled by a continuous functional form for V (x) 
and m(x), which are free from the above mentioned defect of vanishing mass at infinity. A new 
combination of rational forms will be chosen for both functions inside a region, but outside they 
will be considered to be constant, that is to say, our EM Hamiltonian is asymptotically equivalent to 
conventional constant-mass Hamiltonian. Our other purpose is to study the properties of bound states 
of this model and to compare the results with the constant-mass problem. The precise model chosen 
in the present work for m(x) and for the potential is motivated by the following property: it allows 
the transformation to the auxiliary spectral problem of a conventional (constant mass) Schrodinger 
equation, which is solvable. In this sense this variable mass problem can also be characterized as 
solvable. 




(1.2) 
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-x x 
Figure 1. A schematic representation of the model: a variable mass and a potential well in the interval 
(— xq, xo), and outside of this interval both are constant. 

The structure of the paper is as follows. In Sec. 2 we will introduce our model and derive 
the appropriate matching conditions. Applying these conditions the equation for the energies of the 
bound-states will be obtained (the details of some of the issues dealt in this section will be provided in 
the |Appendix A| for readers' convenience). In Sec. 3, we will map the whole problem in an auxiliary 
constant-mass Schrodinger picture and obtain the correct physical range of bound-state energies of 
our EM potential. A limiting case will be considered in Sec. 4, where the potential becomes the 
well-known harmonic oscillator. In Sec. 5, we will solve numerically the transcendental energy 
equation for bound states and examine the spectral properties for a wide range of the parameters. 
Some characteristic differences from the constant-mass case will be noted here. Sec. 6 will contain 
a description of some interesting features of bound-state wave functions, obtained by the analysis of 
the numerical results of the previous section. Finally, we will end with our conclusion in Sec. 7. 



fv{x) = ^3 + - 2 + Ck' 2 x 2 + D. (2.2) 



2. A model Hamiltonian and the energy equation 

2.1. The model 

To start with, we will introduce a hypothetical sample shown in Fig.^ consisting of two materials 
which produces a variation of mass in a stretchable region (— Xq,xq). We are interested in this 
exposition to study the bound states of a particle in the presence of a local potential V(x) 

ffv(x), \x\ < Xq, 

v h (2.1) 
fv{xo) = Vq = constant, \x\ > x 

where the function fv(%) is chosen to be 

Ak 2 Bk 2 
1 + x 2 1 + k' 2 x 2 

The potential function depends on three classes of parameters : i) the ordering parameters 
a, /3, 7 € K of the kinetic term 11.21 . such that a + f3 + 7 = — 1 , ii) the elliptic modulus parameter 
k 2 6 (0, 1) or the complementary modulus k' 2 = 1 — k 2 , and iii) the "Lame parameters" /j,i/6N 
(this terminology will be explained shortly). The four constants A, B, C, D are expressed in terms 
of these parameters as 

A = 4(l + /3- ?? )-( M +i) 2 , B= ( I /+i) 2 -4(l + / 3-7 ? ), 

C* = 2[8?7-5/3-6] -2, D = (1 + /3)(3fc 2 + 2) - 4rjk 2 + ^ - 1. 

In the above expressions the quantity r\ stands for 77 — 1 + (3 + a(a + 13+1). The point xq depends 
on k by xo = 1/Vk* and the mass function is given by 

. [ fm{x), \x\ < X , 

m(x) = < (2.4) 

{ fm(xo) = m = constant, \x\ > x 

where 

f m (x) = [(l+z 2 )(l + fc'V)]~ 1 . (2.5) 
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It should be stressed that the choice of this particular model is made because it not only provides 
a continuous mass-function with physically reasonable non-vanishing limit at infinity, but also leads 
to well-known associated Lame equation, as we shall show subsequently. Thus, our task now is to 
solve the eigenvalue equation II. 1> for the general kinetic energy operator II. 2\ with the potential 
V(x) and mass m(x) given by I2.1l i-( l2~5l . Choosing for our convenience the scale K 2 = 2, this 
equation may be written as 



VAx) + 



m{E-V(x)} 



1 + (3 m" 



V 



ip(x) = 0, 



(2.6) 



where throughout this article prime and dot will denote derivatives with respect to x and z, 
respectively. 



2.2. Matching conditions 

To derive the correct matching and boundary conditions for the envelope wave function i/j(x), we 
notice that the presence of m"/m in 12. 6> produces a ^-discontinuity at the two junctions x — ±Xq. 
Thus, ip' must have a definite jump at these junctions to balance these singularities in the EM 
eigenvalue equation 12. 61 . To calculate precisely these jumps, we will express m(x) and V(x) in 
terms of the Heaviside O-function 





f 1 ' 


x > 


Q(x) = < 


1/2, 


x = 




, o, 


x < 



[G'(x)=6(x)] 



(2.7) 



as 



m(x) = to E(x) + f m (x) [1 - E(x)} . 



where 



V(x) = V E(x)+f v (x)[l-E(x)}, (2.8) 

E(x) = Q(-x+)+Q(x-), x±=x±x . (2.9) 

It may be noticed at once that, due to cancellation effect, m! will not contain ^-discontinuity. 
Thus, we get 

m'(x) = fUx)[l-E(x)}, m"{x) = C{x)[l-E{x)]-f' m {x)E'{x). (2.10) 



It is obvious that the first order derivative term in (12. 6> could be eliminated by introducing the 
transformed function 



c/)(x) = ip(x)/ y / m(x) . 
The EM eigenvalue equation (12. 6> reduces to 

4>"(x) 



c , , , dm" 
m{V(x) -E} + t- 

2 TO 



<t>'(x) 



x +e 



fm{fv - E} + TT-jT 1 



T- 7 ? — 



, , fm 

7 - V -T- 



<j>{x). 

Integration of equation ( 12. 12t with respect to x in the interval [xq — e, xo + e] yields 

4>{x) 



(2.11) 



(2.12) 



+ m (V -E) 



xo+e 



P 



xn + t 



f 

■I III 



x)dx-'- / J -f^5{x-)ct>{x)dx. 

2 Jx -e Jm 



(2.13) 



Let us notice that for our model the wave function ip(x) must be continuous^ on the full line, for 
otherwise ip" (x) would have stronger singularity than that from the term to" in equation 12. bt . It then 

% The case of discontinuity in both V (x) and m(x) (e.g. potential-mass step [ 19 1 or a quantum well with mass mismatch) 
corresponds to discontinuity of M){x) 1511 . 
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follows from \2.\ It that <fi(x) is a continuous function. Hence all integrands in fl2. 1 3i are continuous 
except the last term. Thus letting e — > 0, we will obtain the jump at x — x : 

P /mOo) 



Ad>'\ 



-Mxq). 



2 fm{x ) 

Proceeding similarly for the other junction x = —Xq\ 

_ Pfmi-X*) 



(2.14) 



-6{-xo). 



E » 2 f m (-x ) 

We thus derive the following conditions for the envelope wave function if> and its derivative 
a) ip( x ) is continuous , 



(2.15) 



x— ±a?o 

and for the bound state wave functions 



2 f m (±x ) y/m(±xo) ' 

/oo 
\ij) n (x)\ 2 dx — 1, (normalizability) . 
-oo 

2.3. Expressions for wave functions 

The EM eigenvalue equation (12. 12i takes the following form in the three regions 

K 2 (f>(x), \x\ > Xq 



(2.16) 



(2.17) 



4>"(x) 



fm(x){f V (x)-E} + §^ 



(X), \x\ < Xq, 



in which we have used the abbreviation 

2 -m (V -E). 



(2.18) 



(2.19) 



In this exposition we are primarily interested in bound states and thus the quantity k in ( 12.191 1 
will be always real and nonzero. The scattering states for which E >Vq will be discussed at length 
elsewhere. The acceptable solutions of 12. 181 for \x\ > xq are 



(2.20) 



JVTe , X < -x Q 

<j>{x) = I 

[ N + e~ KX , X > Xq, 

where k is taken as positive square root of ( 12. 19t and the constants W 1 * 1 have to be determined from 
the conditions a)-c). In the intermediate region the reformulation of equation ( 12. 1 8i is useful to 
solve the spectral problem. The crucial observation is that this equation may be transformed into the 
well-known associated Lame equation [52-57] 



X(z) - k 2 



fi(p + l)sn 2 z + v{v + 1) 



dn 2 z 



E 

fc2 



(2.21) 



by means of the following changes of variables 

x(z) = snz/cnz, \( z ) = [/m(^(^))] 1/4 <P(x(z)) . (2.22) 

In the above equations snz = sn(z, fc), cnz = cn(z, k) and dnz = dn(z, k) are three Jacobian elliptic 
functions of real modulus k. Since the parameters /i, v enter in the associated Lame potential 12.211 . 
we will reserve the terminology "Lame parameters" for them. Note that the equation J2.2H is of 
period K or 2K according as fi = v or fi ^ v, where K(k) = j^ 2 dr/y/l — k 2 sin 2 r is the 
complete elliptic integral of second kind. The Lame parameters can be chosen as any integral pair, 
but for simplicity we will choose /i = v = 1. In this context it may be emphasized that the associated 
Lame equation (12.2 1 i for /i = v can be mapped via coordinate transformation z = (1 + k')z 
to ordinary Lame equation with a different modulus parameter k = (1 — fc')/(l + k') and the 
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energy variable E = — //(/i + l)fc + E/(l + k') 2 through the use of the relation sn(z, k) = 
(1 + fc')sn(z, k)cn(z, k)/dn(z, k). However in our auxiliary variable z, equation ( 12. 2\\ represents 
associated Lame equation, and for /i = v it becomes if -periodic. This means that we are going to 
consider this equation in a single period (—K/2, K/2). 

At this stage it is worth mentioning that the region (— xq, xq) of material 2 must be so chosen 
that it may be stretched as large as we please, but it will be always finite (it must be surrounded by 
material 1). On the other hand, it can not be shrunk to a point due to the presence of material 2 (see 
Fig.Q. For this reason we need to exclude the points z = ±K from the domain of equation J2.2H . 
as it can be easily verified from the transformation (I2.22> that these points lead to x = ±oo. For 
definiteness we have chosen the interval z € (—K/2, K/2) which just corresponds to x 6 (— xo, xo) 
where xq = 1/ \fk/. Let us remark about two well-known limits k — * 1 and 0, which are usually 
considered for elliptic functions. Note that as k — > 1, the complementary modulus k' — ► and 
consequently xq goes to infinity. Thus k — ► 1 limit is not allowed in our model. But the other 
limit k — > may be considered as it allows the shrinking of the region (—xq,xq) up to a finite 
interval (—1,1). The general solutions of equation ( 12. 211 for arbitrary energy E, which we need, 
was obtained only recently in Ref 1 571 1581 . Here we will not describe the method of obtaining these 
solutions, but for readers' convenience, we have included a self-contained brief introduction about 
elliptic functions in |Appendix A| The two linearly independent solutions + of ( 12. 2\\ are (for three 
exceptional cases see below) 

*. a w = - p - iK> ^ * * ^ • (2 - 23) 

where ai{E) are to be determined from the equation p(ai) = Ci, c% being zeros of the following 
quadratic equation (see Eqs. (4)-(5) in Ref |57|) 

c 2 + (E - 4 + e x )c + (3ei - e 2 - 2 ei e 3 - e x E) = . (2.24) 

Here, e, are always real (e\ > e 2 > 63) defined by p{u>j) — e,; u>\ and uj 3 are half-periods of 
Weierstrass elliptic function p(z) (see |Appendix A| for more details). For the numerical convenience, 
here we have chosen the scale e\ — e 3 = 1 so that u>i = K, 003 = iK', L02 = u)\ + UJ3, K'(k) = 
K{k'). Thus, the solution in the intermediate region takes the form 

xO) = diXi(*0 + d 2X2{z) ■ (2.25) 

The constants dt, c?2 have to be fixed from the conditions a)-c). Before proceeding to do that, we 
would like to point out that for special values of E = E^\j = 0, 1, 2 both solutions Xi( z ) an d 
X2( z ) become identical and coincide to periodic (or anti -periodic) band-edge solutions X ■ The 
explicit expressions for these solutions 1 53 1 and corresponding values of a, in ( I2.24> are 

X (0) = dnz + fc'/dnz, E^=2 + k 2 -2k', 01 = -a 2 = wi/2 , 

= dnz - k'/dnz, = 2 + k 2 + 2k', 01 = -a 2 = -u 3 + Wi/2 , (2.26) 

^-( 2 ) = snzcnz/dnz, E^ =4, a% = 0J3, a 2 = oj 2 . 

Thus, for E = E^' the two linearly independent solutions of associated Lame equation ( I2.2H may 
be given as 

/ z dr 

It should be emphasized that in our model the situation is different from that of periodic 
associated Lame equation 1 53 , 33] because the equation ( 12. 2H is considered in only one period 
(—K/2, K/2). We will show later that the discrete energy levels of our EM Hamiltonian lie inside 
the allowed bands of constant-mass periodic associated Lame Hamiltonian, and corresponding wave 
functions ip(x) are obtained from (I2.20> and J2.23I . 

+ The notations for the Lame parameters fi, v in Ref. 57 58 are m, i. 



A New Effective Mass Hamiltonian and Associated Lame Equation: Bound States 



7 



2.4. Energy equation for bound states 



The continuity condition and slope-discontinuity requirement J2.161 can be expressed in a single pair 
of equations 



xt 



d 2 xf 



K)e 



—k/ yW 



E Xi ) + d 2 (TmoX2 2~ TO o 



l xf) 



where we have used the abbreviations 

Xf = Xi(±K/2), X f = Xi{±K/2), i = 1, 2. 
Eliminating 3V ± from J2~28l -l l2~29l . we obtain a homogeneous linear system for d\ , d 2 



(2.28) 
(2.29) 

(2.30) 



d x [2? 



1/4 



^xt 



d 2 [2i 



1/4 



-1/4 



4^ 



where the quantity 23 reads 

3 = (2(3+ 1> 



lQ 1/4 03 X r] + d 2 [2my 4 X2 - - m '"Sxj ] = 0, 



-1/4, 



2k. 



(2.31) 



Demanding for non-trivial solutions of d\ , d 2 from the previous system of equations, we have 
obtained our energy equation in the following form 

4m (T 2 - - T+) + 2V^'B(Tf + T+) + S 2 (T+ - Tf) = 0, 



(2.32) 



where 



T? 



xtxf; 



t 



xtx 2 



xtx 2 \ 



T, 1 



xtx 2 - 



1 — Al A2 > ^1 -A1A2 Al A2 ' J 2 — Al A2 ■ (2.33) 

In particular for E ^ E^\j = 0, 1,2, the energy equation \232l may be further simplified by 
inserting the expressions for xt fr° m equation ( I2.23> (see |Appendix A| > and then the energy equation 
reduces to the form 



[Am A 2 _ - 4£- v /ro^4_ + B 2 ] T i - [4m A 2 + + A r B^ / m^A + + S 2 ] T+ = 0, 



where A± are given by 



A+ = ±- 



' ' fp(f + m 3 ) Tp(ai) p(f + ^ 3 ) 



p(<H 



• w 3 1 - ei 



(2.34) 



(2.35) 



Hence, similar to the well-known textbook example of the finite square well potential, we have 
to deal with the transcendental energy equations J2.32> and J2.34l >. It is clear from the structure of 
equation J2.34l > that for T\~ = some roots will be given by the equation A- — A + = 0. But one 
may notice from J2.35l > that this equation will be true if and only if ^2 p( a i) = 0. The latter condition 
is, however satisfied only for band-edge energies E = E^\ mentioned in the equation \2.26\ . This 
clearly implies that the roots E = of (I2.34> have to be discarded, because for these values the 
correct energy equation is J2.32> . where Xi & are given by J2.26> and J2.27i . It may be mentioned 
that we have checked numerically for different values of the potential parameters that the energy 
equation J2.32i is not satisfied for E = E^> . This means that for those cases the band-edge energies 
of periodic associated Lame Hamiltonian are not roots of the energy equation for bound states of our 
EM Hamiltonian. Thus, we will solve numerically the equation J2.34i for E, and the roots, if they 
exist, may lie in principle inside the allowed and forbidden bands. However, the numerical results 
could be well-handled if we have, a priori, some analytical insight from the theory. To achieve this 
insight we have to fall back upon the auxiliary constant-mass Schrodinger equation of which we 
know very well the structure of the spectrum. In the next section, we will obtain an image of our EM 
model in this conventional Schrodinger language. 
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3. Auxiliary constant-mass equation 

Our purpose is to extend 12. 211 onto the whole z-axis in the form 

-x(z) + V(z) X (z) = E X (z), (3.1) 

which can be viewed as a constant-mass Schrodinger equation for a single particle of unit mass 
(according to our chosen scale ti 2 = 2). Of course, equation 13. 1> plays the role of auxiliary equation, 
and our aim is to extract general information about its spectral properties, which could be used as a 
guide in our numerical procedure for physical model II. 11 1. The observation here is that this can be 
achieved by extending the transformations I2.22> to the entire z-axis 

\/m(r)dr, xO) = [m(x(z))}~ 1/4 ip(x(z)) . (3.2) 

But there are two subtle points which should be taken into account to obtain the correct expression 
for the Schrodinger potential V(z). In the first place, since m(x) is continuous, the new coordinate 
z(x) should also be continuous. It will then follow from 13. 2> that the constant-mass wave function 
X(z) is also continuous function of z. One can achieve this by exploiting the arbitrariness in the 
indefinite integral in (13.21 . Indeed the explicit relation between z and x is 

-oo < z < -K/2, 

-K/2 < z < K/2, (3.3) 
K/2 < z < oo, 

where A± = ±(K/2— y/moXo). The second point is even more fundamental concerning the nature of 
m! and to" in the whole region iei We have mentioned in Subsec. l2~2l that due to the cancellation 
effect to' will contain 0-discontinuity and consequently to" will produce ^-discontinuity at the two 
junctions x = ±xq. Noting that S(bx) = 6(x)/\b\, the Schrodinger potential V(z) may be expressed 
in the following form 

V(z) = V S(z) + f v (z)[l - S(z)] + (fl + (1 + k') [S(z+) + 6(z-)] , (3.4) 




where 

fv(z) = 2k 2 



sn 2 z + 



dn 2 z 



E(z) = e(-z+) + Q(z-), z± = z±K/2. (3.5) 



The point to be noticed here is the /5-dependence of the shape of V(z) at the two junctions 
z = ±K/2. It should be kept in mind that although the auxiliary constant-mass Schrodinger wave 
function X (z) is continuous, its derivative is not. In this sense the auxiliary Schrodinger equation 
13. 1> is not exactly in conventional form. One may check readily from the matching condition I2.16> 
that 

A [*(*)] _ =(p+^) (1 + A;')X ± - (3-6) 



z=±K/2 \ 2 

To understand the correct range of E for bound states, we will consider separately three cases: 
> -1/2, [3 = -1/2 and /3 < -1/2. 

• Casel)/3 > -1/2 

In this case there will be a well, determined by one-period associated Lame potential fv(z) in 
the region (—K/2, K/2), bounded by two i5-barriers, and outside this region a constant potential 
Vq. Three possible situations may arise (see Fig.|2jl. Clearly there may be bound states in the 
range (f v )min < E < V if and only if (fv)min < V . 

• Case 2) (3 = -1/2 

In this case the <5-barrier will disappear and for bound states the conclusion is the same as before. 
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8 barriers 




-K/2 K/1 2 

Figure 2. Three possible positions of the constant-mass Schrodinger potential V(z) are shown for 
/3 > —1/2. For /3 = —1/2 5-barriers disappear, and for j3 < —1/2 there are two 5-wells at the 
junctions. 



• Case 3) (3 < -1/2 

In this situation the shape of the associated Lame well provides no restrictions from below on 
the bound state energies due to the presence of two <5-wells at the two junctions z — ±K/2. We 
thus conclude that the bound states may exist in the range — oo < E < Vq. 

In the next section we will consider the limit k — > and will write the explicit energy equation 
for this limit. Before concluding this section it may be mentioned that in the k — > limit, fv{z) = 
and V(z) reduces to 

V(z)->V E(z) + (2p + l)[6{z+)+6(z-)] , (3.7) 

where z± — > z ± tt/4. Thus, for (3 > —1/2, (Case 1 and 2 above) bound states may exist in the range 
< E < Vo and for (3 < -1/2 the possible range is -oo < E <V Q . 



4. A limiting case : k — > 

It is well-known that the Jacobian elliptic functions degenerate into trigonometric and hyperbolic 
functions in k — > and k — > 1 limits. We have akeady pointed out that k — * 1 limit is prohibited* in 
our model, as this limit corresponds to the nonphysical situation of vanishing mass for large |x|. In 
the k — > limit, the intermediate region shrinks up to the interval (—1,1) inside which the harmonic 
oscillator well V(x) and the mass function m(x) are 

T/ , . f Cx 2 + D, \x\ < 1 . . f (1 + x 2 )- 2 , |x|<l 



V , \x\ > 1 ' " w \ 1/4, |x| > 1 

In the two infinite regions — oo < x < —1 and 1 < x < oo, the solutions will be given by ( 12.11 
and J2.20I . where k, vtlq and Vq take the limiting forms 



K — * t/Vo-E/2, Tno-fl/4, F 5 + 8 [/3 + 2a(a + /3 + 1)] . (4.2) 
In the intermediate region, we have obtained the free-particle equation 

-x(z) = E X (z), «e(-7r/4,7r/4), (4.3) 
where the transformations J2.22i now reduces to 

a; = tanz, ^(z) = (sec z)i/j(x(z)) . (4.4) 

* But this limit becomes interesting if different boundary conditions are adopted 59 1. 
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K-e KX /2, x < -1, 

^givStan (as) _|_ tan~ 1 (a;) 
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X < 1, 



(4.5) 



VTTx 2 

Ji+e- KX /2, at>l. 
As before, imposing the matching conditions, we obtain homogeneous linear systems for d\ , d 2 

die-^/^B - iVE) + d 2 e l ^l A [¥> + iVB] = . 
Thus the energy equation for bound states (E < Vq) in the k — > limit is as follows 

(CB 2 - B)sin(|"/B) + 2SV!Ecos(|\/s) = 0, (4.6) 

where 23 is given by 12.3 It for k' — > 1. One may notice at once that _E = is a trivial root of equation 
14. 6> . but it must be rejected, because the solution ^(s) for the intermediate region (—1, 1), given in 
14. 5> . is valid for £7 7^ 0. Indeed, for E = one may rewrite the solution as 

V>(x) = [d\ tan -1 x + d 2 ], — 1 < a; < 1, (4.7) 

the expressions in two semi-infinite regions being the same as in J4-.5I . In this case the matching 
conditions give the following constraints for the existence of zero energy root 

A 1 2 



V = (2/3+ l) 2 , /3 < 



or V = 



2/3 + 1 + 



1 



(4.8) 



and correspondingly in ( 14.71 either cfi = or ^2 = 0. We have checked that the constraints 14.81 are 
not satisfied for different cases, so that zero energy state does not exist for them. 

5. Numerical results for bound states 



In this section we will solve numerically the energy equations ( 12.341 [and (I4.6i for k — > limit] in the 
range V m i„ < E < Vq, where V m in denotes the minimum of auxiliary potential V(z) given by d3.4i 
[and ( 13. 7H . It may be mentioned that the energy equations depend on two classes of parameters : i) 
ordering parameters a, (3, 7 connected by a + (3 + 7 = —1, and ii) elliptic modulus k 2 (0 < k 2 < 1) 
or the complementary modulus k' 2 = 1 — k 2 . Thus we have examined the roots of the energy 
equations for bound states as a function of both these two classes of parameters. Our strategy is to 
vary these parameters in such a way that it will cover some of the special forms of kinetic energy 
operator II. 2> mentioned in the introduction. 



5.1. One parameter family of kinetic energy operator for a = 7 
In this case kinetic energy operator II. 2> reduces to 

Tem(x) = X - (m a pm pm~<) , 2a + f3 = -l. 
We will first obtain E as a function of k 2 £ [0, 1) for the two values of (3: (3 = — 1 and (3 = 0. 



(5.1) 



Case/3 = -1 

For the choice /3 = — 1, the kinetic energy operator 15. 1> takes the following form 

Tem(x) = - (p—p 
2 \ m 



(5.2) 
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Table 1. The bound state energies are calculated for /3 £ [—2, 2] and a = 7, fc 2 = along with 
M), Kni n (xq = 1). The bound states are observed only in the range [—0.4, —0.1]. 







V 


E 


-0.10 


0.80 


0.96 


0.94 


-0.20 


0.60 


0.84 


0.81 


-0.30 


0.40 


0.64 


0.62 


-0.40 


0.20 


0.36 


0.35 



Although in this case the upper and lower limits for bound states are 

9fc 2 (l-fc') , k 2 

and V m i n = —00, since j3 < —1/2 [see equation J3.4H . no bound states are obtained in this range. 
Case (3 = 

This choice of (3 corresponds to the kinetic energy operator 

Tem(x) = I l-^p 2 -^] . (5.4) 



2 \ 0n yrn , 
Here, the range is 

V °- 5 -WTW + 1 + ^ V mm = 2P. (5.5) 

In the interval (V m i n , Vo) there exist no bound states for k 2 € [0, 1). 

Next we will study the bound-state energies as a function of (3 € [—2, 2] for the two cases: 
k 2 = and k 2 = 0.5. 

Case k 2 = 

The upper limit Vo has a parabolic dependence on (3, given by 

V = 1 - 4/? 2 , (5.6) 
and the lower limit is 

V mm = { I (5.7) 

[ -00, -2<(3< -|. 

Our numerical calculation shows the existence of one bound state near the top of the well for 
-0.4 < < -0.1 (see Table 1). 

Case k 2 = 0.5 

In this case also the upper limit Vo has a parabolic dependence on j3, given by 

while the lower limit is given by 

f 1, ~\<P<% 
V mm = { I (5.9) 

[ -00, -2 < (3 < -\. 

No bound states exist in this range. 
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Table 2. Bound state energies are provided for k 2 E [0, 1) where a = —l,/3 = 7 = 0. The point rrrj 
and Vo, V m i„ for V(x) are shown in each case. Bound states start to appear from fc 2 < 0.90. 



k 2 


x 


'min 


Vo 


E 


0.90 


1.78 


0.32 


2.64 


2.63 


0.80 


1.50 


0.40 


2.88 


2.59 


0.70 


1.35 


0.48 


3.13 


2.51 


0.60 


1.26 


0.55 


3.38 


2.43 


0.50 


1.19 


0.62 


3.65 


2.35 


0.40 


1.14 


0.70 


3.91 


2.28 


0.30 


1.09 


0.78 


4.18 


2.20 


0.20 


1.06 


0.85 


4.45 


2.12 


0.10 


1.03 


0.92 


4.73 


2.04 


0.00 


1.00 


1.00 


5.00 


1.96 



5.2. Two parameter family of kinetic energy operator for a ^ 7 

We will find E as a function of k 2 for a = — 1, j3 = which yields following kinetic energy operator 



Tem(x) 
The bounds are given by 



1 



1 



-P +P 



V a = 



3fc 2 (l-2fc') 



4/c' + 1. 



V m 



= 2k 1 



(5.10) 



(5.11) 



2(1 + */) 

We have obtained one bound state of the potential (see Table 2) for < k 2 < 0.90. The effective 
mass potential V(x) maintains the shape of the well for all k, and in particular for k — it becomes 
well-known harmonic oscillator inside the interval (— xq, Xq), the minima being at x = 1, It is 
interesting to observe that the bound state exists inside the first allowed band of periodic constant- 
mass associated Lame potential (see Fig.|3j- 




Figure 3. Bound state energies as a function of k 2 are shown by dots for the case a = — 1, 
j3 = 7 = 0. The lower and upper gray curves represent first two band-edge energies 
E^',j = 0, 1 for the corresponding constant-mass periodic associated Lame potential. 
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The dependence of bound-state energy E on f3 will now be calculated again for two cases: 
k 2 = and k 2 = 0.35, while (3 will vary in the same interval [—2, 2] and a will be fixed as a = —1. 

Case k 2 = 0,a = -1,7 = -/3 

Here the upper limit Vo is linearly dependent on (3 and is given by 

V Q = 5-8/3, (5.12) 
and the lower limit V mm is same as given by d5.7i . The results, shown in Fig.0] are the following 




-2 -1.5 -1 -0.5 0.5 1 

Figure 4. Bound state energies (dots) as a function of j3 for a ^ 7, a = — 1, k 2 — 0. 



• No bound states for 0.5 < (3 < 2, 

• One bound state for -0.1 < (3 < 0.4, 

• Two bound states for -1.0 < (3 < -0.2, 

• Three bound states for —2.0 < (3 < —1.1. 

One can notice that for (3 = 0.4, V(x) = 1.8 is a constant potential (see Table 3), but still there 
exists one bound state. However, we have already explained in Sec.[3]that this apparently strange 
behavior is perfectly consistent with constant-mass Schrodinger picture. In this context it may be 
mentioned that the existence of bound states for constant potential was also noticed in Ref. |38| for 

m{x) = sech 2 ga;. 



Case k 2 = 0.35, a = -1,7 = -/3 

This case also corresponds to linear dependence of Vo on j3 
Vo = 4.05-5.16/3, 
and the lower limit is given by 

' 0.7, /3>-l/2, 
-00, /3<-l/2. 
The results shown in Table 4 are : 

• No bound states for 0.5 < < 2, 

• One bound state for -0.2 < (3 < 0.4, 

• Two bound states for —1.2 < (3 < —0.3, 



Vmin — 



(5.13) 



(5.14) 
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Table 3. The bound state energies as a function of j3 G [-2, 2] for a = -1 ^ 7, k 2 = are given 
with V m i„ and Vo (xo = !)■ The asterisk in /3 indicates that is constant potential. 



/3 


^rmn 


Vo 


E 


0.40* 


1.80 


1.80 


1.69 


0.20 


1.40 


3.40 


1.91 


0.00 


1.00 


5.00 


1.96 


-0.20 


0.60 


6.60 


1.96; 6.36 


-0.40 


-0.20 


8.20 


1.94; 6.83 


-0.60 


-0.20 


9.80 


1.91; 7.00 


-0.80 


-0.60 


11.40 


1.86; 7.03 


-1.00 


-1.00 


13.00 


1.79; 6.97 


-1.20 


-1.40 


14.60 


1.72; 6.86; 14.14 


-1.40 


-1.80 


16.20 


1.64; 6.70; 14.61 


-1.60 


-2.20 


17.80 


1.54; 6.48; 14.79 


-1.80 


-2.60 


19.40 


1.43; 6.24; 14.81 


-2.00 


-3.00 


21.00 


1.31; 5.95; 14.72 



Table 4. The bound state energies as a function of /3 € [—2, 2] for a = — 1 7^ 7, k 2 = 0.35 are given 
with Vmi n and Vo (xo = 1.11). 







Vo 


E 


0.40 


1.96 


1.98 


1.98 


0.20 


1.35 


3.01 


2.20 


0.00 


0.74 


4.05 


2.24 


-0.20 


0.13 


5.08 


2.23 


-0.40 


-0.48 


6.11 


2.20; 5.67 


-0.60 


-1.09 


7.14 


2.16; 5.86 


-0.80 


-1.70 


8.18 


2.10; 5.89 


-1.00 


-2.31 


9.21 


2.04; 5.83 


-1.20 


-2.92 


10.24 


1.96; 5.71 


-1.40 


-3.53 


11.27 


1.87; 5.54; 11.04 


-1.60 


-4.14 


12.31 


1.77; 5.34; 11.35 


-1.80 


-4.75 


13.34 


1.66; 4.96; 11.43 


-2.00 


-5.36 


14.37 


1.53; 4.82; 11.40 



• Three bound states for —2 < j3 < —1.3. 

Once again we have observed that the ground state lies inside the allowed band and higher 
excited states are inside the continuum (see Fig.|5j for the constant-mass periodic associated Lame 
potential characterized by /i = v = 1; k 2 = 0.35. 

6. Some features of bound-state wave functions 

In this section we will present some interesting properties of the wave functions i/j(x) that can be 
derived from the numerical results of the previous section with very good accuracy. A part of these 
properties, which seems to be difficult to prove analytically, does not concern the specific physical 
problem. We recall that ip(x) is given by (up to the normalization factor) 

!m Q ~^ 4 x~ e K ( Xo+x \ —00 < x < —xo , 
rt 4 (x)x(z(x)), -x <x<x , (6.1) 
m ~^ 4 x + e K ( x °~ x \ xo < x < 00 , 



A New Effective Mass Hamiltonian and Associated Lame Equation: Bound States 

12rE 



15 




-1.5 -1 



0.5 



Figure 5. Bound state energies (dots) as a function of f3 are plotted for a ^= 7, a = — 1, 
k 2 = 0.35. The gray lines represent band-edge energies of corresponding associated Lame 
potential. 



where x( z ( x )) = X2 + (di/d 2 )xi(z(x)) and x ± = x(±^/2). Note that for k ^ 0, 
Xi(z(x)) are given by J2.23l > for x{z) = snz/cnz, and for k — 0, Xi,2(2 : (a;)) = exp[±iV A E2;] for 
= tan z. The quantity k is to be computed from J2. 19i for E = E n , E n being the bound-state 
energies of our EM Hamiltonian. In the following we will use the notation 



«n = V m o( v o - E n ), n = 0, 1, 2 . 



(6.2) 



6.1. Ground state 

The two quantities <Zi and a 2 in ( 12.231 1 are complex and related by 

<ii = —a^ . (6.3) 

The real and imaginary part of Xi (z) are respectively odd and even functions. Moreover Xi ( z ) satisfy 
following interesting relations with very high accuracy 

d 2 



xiO) 



dx 



X2(-z), Xi(z) = -X*(- Z )i 



(6.4) 



where the ratio d-z/dx is a positive number. It will then follow that x( z ) i s an even function and 

X(z) = x(-z) cx lm[xi(z)} ee v{z) . (6.5) 
The ground state wave function is 



e ml /4 v( f ) e «o(*o+z) j 
M*) = { QofU\x)v{z{x)), 

Q ml /4 v (f ) e K °( x °- x \ 



X < —Xq 

\x\ < x , 

X > Xq . 



(6.6) 



It is straightforward to compute the normalization constant Co from (12.171 



1 



v 2 (K/2) 

y/Vo — Eq 



V fm(x)v 2 (z(x))dx 



1/2 



(6.7) 
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6.2. First excited state 

The quantity ai is purely imaginary while the other quantity a 2 is complex and they are related by 



Im[ai] = Im[a2] ■ 
The relation between X\(z) and \2 (z) are 



xiO) 



'di 



where d 2 /d\ is a positive number. Clearly this implies that 

x{—z) = — x( z ) = complex . 
We have checked numerically with the same high accuracy that 

Re[x(*0] oc Im[x(^)] , 
and of course no degeneracy exists for E = E\ . The wave function is 



6.3. Second excited state 



-eim; /4 Re[x + ]e Sl ( I »+ I l, x < -x , 
e 1 f 1 J 4 (x)R S [ X (z(x))}, \x\<x , 
eiTOo /4 Re[x + ]e Kl ( a:o - a: \ x > x a . 



(6.8) 

(6.9) 

(6.10) 
(6.11) 

(6.12) 



The two quantities a\ and a 2 have same properties as in the above case. The relation between xi{ z ) 
and X^{z) reads in this case 



Xi{z) = j^X2{-z), 



Thus, here also x(z) is complex and 

x(-z) = x{z), 

but 

Re[x(z)} cx lm[ X (z)} . 
The wave function is 

e 2 TOo /4 Re[x+]e K2 ( ;ro + :E \ x < -x , 
Mb) = \ Q 2 fU\x)Re[ X {z{x)% \x\ < x , 
k e 2 m 1 Q /4 Re[x + }e K ^ X0 - x \ x > x . 
The normalization constants C„, n = 1, 2 may be expressed as 



(6.13) 

(6.14) 
(6.15) 

(6.16) 



(Re[x + ]) 2 



^UJx){Re[x{z(x))]fdx 



1/2 



(6.17) 
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6.4. The limiting case for k — > 

In the limiting case for k — > the relation between x and z is very simple [see equation J4.4H . and it 
is convenient to express the corresponding properties in the variable x. Here we will denote Xi(z(x)) 
and x(z(x)) by Xi{ x ) an d x( x )- The relation between xi( x ) an d X2( x ) is 

Xi(-x)= X 2(x). (6.18) 

Let us introduce the parameters 

O n = j\ZK, n = 0,1,2. (6.19) 

Ground and second excited states (k = 0): 

In both cases the ratio efe/di is unity. It then follows from the definition of x( x ) that it is even 
function given by 

x(-x) = x(x) = 2 cos(\/£? tan -1 x) . (6.20) 
The wave functions are 

e n V2cos9 n e K " (1+x \ x<-l 
^ n (x) = { -^^(V^tan- 1 ^), |z| < 1 (6.21) 



n/1H 

e„%/2cos0„e K "( 1 - ;r ) x<-l 



for n = 0, 2. 



First excited state (k — 0): 

The ratio dz/di is equal to —1 so that x( x ) is an odd function given by 

x (x) = -x(-x) = -2isin(V£'tan" 1 x) . (6.22) 
The wave function is 

{-Ci\/2siri6>ie Kl(1+:E) , x < -1 
^f|p ^(VT^tan- 1 ^), \x\ < 1 (6.23) 
eiV2sin6»ie Kl ( 1 - a; ), x<-l, 
In all the three cases the normalization constants C„ may be expressed as 

' 1 + (-!)" cos2fl w | ^ n sin 20, 



1 



1/2 

(6.24) 



7. Conclusion 



In this article we have proposed a new solvable model wherein the potential and effective mass are 
rational functions of spatial coordinate. The novel feature of our model is that it may be mapped 
to well-known periodic associated Lame potential in constant-mass scenario. This fact is clearly 
responsible for solvability of the model. It may be mentioned that the list of solvable models in 
effective mass framework is rather short. Our work has definitely enhanced this set by introducing a 
wide class of rational potential and mass functions. The important difference of our work compared 
to recent efforts is that we consider the variation of mass inside a finite region, and both potential and 
mass are constant outside. The advantage of considering variation of the mass inside a finite interval 
instead of full line is that the mass remains finite and non-zero everywhere, as it should be. 

We have examined the bound-state spectrum for different values of the ordering parameters 
a,f3,j and elliptic modulus fc 2 , and have discussed the properties of the corresponding wave 
functions. It is observed for both a — 7 and 0^7 cases that the ordering parameter (3 has a critical 
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value above which no bound state exists. But number of levels increases with decreasing values of 
(3. Some peculiarities have been noted in the spectral properties [e.g. the existence of bound states 
for constant potential V(x)] in contrast to the conventional situation where mass is constant. The 
qualitative observation is that the discrete energy levels for bound states of our EM Hamiltonian for 
nonzero k lie only inside the allowed band and in the continuum for corresponding constant-mass 
periodic associated Lame Hamiltonian. Moreover, as a by-product of our numerical procedure we 
have found some new curious relations for the solutions Xi( z ) °f associated Lame equation, which 
seem to be universal. For example the form of the second relation of (16. 4> for the solutions Xi{ z ) at 
E = Eq (inside the allowed band of corresponding associated Lame potential) is not related with our 
physical problem. 

Some direct generalizations of our model are possible. For instance, the application of SUSY 
transformations to enlarge the class of rational models and the inclusion of scattering states with 
E > Vq are under consideration. Further it will be interesting to investigate the validity of the above 
mentioned relations between the solutions of associated Lame equation for other values of E. The 
special exploration of k — * 1 limit of our model also seems to be promising 1 59 1. 
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Appendix A. 

In the following we will mention basic definitions and the relations involving elliptic functions, which 
we use in the text (for more details, see [59-61]. Consider two real numbers k 2 and k' 2 such that 

fc 2 e(0,i), k' 2 = i-k 2 . (A.l) 

These two numbers are called respectively elliptic modulus and complementary modulus, and are 
basic parameters in the constructions of elliptic functions. The amplitude function is defined by 

dr 

tp(z,k) = am(z,fc), z(cp,k) = - == . (A.2) 

Jo y 1 — kP" sin t 

From here the three Jacobian elliptic functions are defined by 

sn(z, k) = sirup , cn(z, k) = cos <p , dn(z,k) = dcp/dz . (A. 3) 

These functions are called sine-amplitude, cosine-amplitude and delta-amplitude respectively. For 
simplicity, in equation J2.2H in the text, and also in the following, we have suppressed the explicit 
modular dependence and write simply sn z, cnz, dnz. These are doubly-periodic functions of periods 
(4:K, 2iK'), (4K, AiK') and (2K, AiK') respectively, and are usually defined for a complex variable 
z. Nevertheless in our case z is always real. The quarter-periods K and K 1 are the real numbers 
given by 

K(k) = K = z(n/2, k), K'(k) = K' = K(k') . (A.4) 
K is called complete elliptic integral of second kind. Noting the following relations 

cn sn z 

m(z + K) = - — , cn(z + K) = -k'- — , dn(z + K) = - — , (A.5) 
dnz dnz dnz 

sn(z + 2X) = -snz, cn(z + 2K) = -cn z , dn(z + 2 K) = dn z , (A.6) 

we see that the associated Lame potential in equation ( 12. 21i is 2X-periodic or X-periodic, according 
as [a ^ v or /i = v. The function snz is odd with a simple zero at z = 0, while cnz, dnz are even 
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functions; cnz has a simple zero at z = K, but dnz has no zeros for real z. The precise value of the 
point x is obtained as x = 1/ \f~ki from the values 



sn lyl 



1 



Some other relevant relations are 

sn 2 z + cn 2 z = 1 , 



K 
cn | — 



VT+¥' 



(A.7) 



dn^z + k 2 sn 2 z = 1 



In the two limits k — > 1 and fc 
trigonometric functions 



(A.8) 

sn'z = cnzdnz , cn'z = — snzdnz , dn'z = — fc 2 sn z cn z . (A. 9) 

0, the elliptic functions degenerate into hyperbolic and 



snz 



fe->i 



tanh z 
sin z 



cnz 



fe->i 

k-^O 



sechz 

cos z 



dnz 



fe-»i 



sechz 
1 



where the quarter-periods go over 

oo 

tt/2 ' 



K 



fe->i 
k-a 



K' 



tt/2 

oo 



Weierstrass elliptic function p(z; g 2} 53) = p(z) is defined by 



1 



(A. 10) 



(A.ll) 



(A. 12) 



(z — 2muj\ — 2ncj 3 ) 2 (2mwi + 2nw 3 ) 2 _ 

where the symbol means summation over all integral values of m, n except m = n = 0, oj\ and 
being half-periods of p(z). In our case, these are defined through 

Wi = K, u 3 = iK'. (A.13) 

and the invariants g 2 , 53 are given by 

32 = ^(fc 4 -fc 2 + l), 53 = ^(fc 2 -2)(2fc 2 -l)(A : 2 + l). (A.14) 

The above choice corresponds to the case when the discriminant A = g\ — 21 g\ > 0, and so the 
three numbers p{u)i) = ej, i = 1, 2, 3 are always real (ei > e 2 > e 3 ), where W2 = Wi + W3. It may 
be mentioned that ei are the three roots of the equation 

4t 3 - g 2 t - .93 - , (A. 15) 

and, for simplicity, we have chosen the scale e\ — e 3 = 1. Weierstrass elliptic function is an even 
function; its derivative p(z) is an odd elliptic function with the same periods and satisfy following 
identity 

3 

p 2 (z)=4l[[p(z)-e i }. (A.16) 

»=i 

The relations between Weierstrass elliptic function and Jacobian elliptic functions, with our 
choice, are 



ei + 



e 2 + 



dn 2 z 



e 3 + 



1 



(A. 17) 



Weierstrass zeta function C{z\92,9?,) = C( z ) an d sigma function a(z; 52,53) = v{z) are quasi- 
periodic functions defined by 



C( z ) = -P(z) , 



a[z) 
a(z) 



and satisfy the following relations 

C(-z) = -C(z) , ((z + 2ui) = ((z) + 2((aJi) , 

(j(-z) = -<t(z), a(z + 2uji) = -a(z) exp[2((uji)(z + uji)] . 



(A.18) 



(A. 19) 
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The addition formulae for these functions are 



pOi + z 2 ) 



p(zi) - p(z 2 ) 
p{zx) - p(z 2 ) 



p( z i) - pM 



C( Zl +z 2 ) =C(zi)+C(^ 2 ) + 



1 p{zx) ~ P(z 2 ) 



(A.20) 



(A.21) 



(A.22) 



2 p{zi) - p{z 2 ) ' 

We will now derive the energy equation d2.34i . Using the relations ( IA. 1 3i — ( |A. 19i and the 
addition formula JA.2H . and noting that 

p(z + 2u)i) = p{z) , p(-z) = p(z), 

p(z + 2ui) = p(z) , p(-z) = -p(z) , 

it is not very difficult to derive the following relations 

xt = A± X f , xf = -A tX2 , (A.23) 

where Xi( z ) ar, d A± are given by 12.231 and J2.35i . Inserting the expressions for \f from ( IA.23t 
into J2.33b . the energy equation J2.34b will readily follow from J2.32l >. 
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